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Abstract 



Optical absorption in a random one-dimensional lattice in the presence of 
paired correlated disorder is studied. The absorption line shape is evaluated 
by solving the microscopic equations of motion of the Frenkel-exciton problem 
in the lattice. We find that paired correlation causes the occurrence of well- 
defined characteristic lines in the absorption spectra which clearly differ from 
the case of unpaired correlation. The behavior of the absorption lines as a 
function of defect concentration is studied in detail. We also show how exciton 
dynamics can be inferred from experimental data by deriving an analytical 
expression relating the energy and intensity of lines to the model parameters. 

PACS numbers: 71.35-Fz; 36.20.Kd; 78.90.-Ft 
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I. INTRODUCTION 



During the last few years, it is being realized that correlated disorder has profound ef- 
fects in quasiparticle dynamics and produces a variety of complex and new phenomena. The 
notion of correlated disorder implies that certain physical parameters in random systems are 
not completely independent within a given correlation length, thus leading to a competition 
between short-range order effects and the underlying long-range disorder. It is well known, 
for instance, that delocalized electronJlFB and vibrationsi appear in one-dimensional ran- 
dom systems with correlated disorder. Following the above line of research, we have recently 
carried out the time-domain analysis of trapping of Frenkel excitons in one-dimensional sys- 
tems with paired correlated traps, randomly placed in an otherwise perfect lattice.i Results 
were compared with the dynamics of one- dimensional lattices whose excited states are also 
Frenkel excitons, but in which transport is affected by the presence of unpaired traps. I We 
found that pairing of traps leads to a slowdown of the trapping rate due to the occurrence of 
larger segments of the lattice which are free of traps. In view of these results, the question 
arises as to whether the fact that disorder is correlated will affect the optical absorption 
spectrum of the lattice. If so, it would be interesting to find a relationship between the 
characteristic absorption lines and the model parameters in order to compare with exper- 
imental data. In this paper we aim to answer these questions by numerical simulation of 
optical spectra. We believe that the interest of this work is beyond the formal study of 
the exciton dynamics: Indeed, our treatment may help describe the optical properties of 
some polymers with active molecules, provided that interchain interactions become much 
smaller that intrachain interactions so that a one-dimensional model can be applied as a 
first approximation. 



II. MODEL 

The Hamiltonian for the Frenkel-exciton problem can be written in the tight-binding 
form with nearest-neighbour interactions as follows (we use units such that h = 1) 

'H = J2 Vkalttk + T ^ {alttk+i + al^^ak). (1) 

k k 

Here and al are the exciton annhilitation and creation operators, respectively. Vk is the 
transition frequency at site k and T is the nearest-neighbor coupling, which is assumed to 
be constant in the whole lattice. The paired correlated disorder is introduced as follows: 
Transition frequencies can only take on two values, V and V, with the additional constraint 
that V appears only in pairs of neigbouring sites (the correlation length is rouhgly the 
lattice spacing). Such a pair of sites will be called dimer defect. In addition, we define the 
defect concentration c as the ratio between the number of sites with transition frequency V 
and the total number of sites in the lattice. 

The line shape I{E) of an optical-absorption process in which a single exciton is created 
in a lattice with N sites is given by (see Ref. ^ for details) 

I{E) = -^j^ dt e-* sm{Et) Im G,(t)^ , (2) 
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where the factor exp(— at) takes into account the broadening due to the Lorentzian instru- 
mental resolution function of width a. The correlation functions Gy^if) obey the equations 
of motion 

ijG^{t) = Y.ii^fiAi). (3) 

with the initial conditions Gfc(O) = 1. The diagonal elements of the tridiagonal matrix Hjk 
are whereas off-diagonal elements are simply given by T. Once these equations of motion 
are solved, the line shape is found using Eq. (0). 



III. RESULTS AND DISCUSSIONS 

We have numerically solved the equation of motion for chains of = 2 000 sites 
using an implicit integration scheme. In order to minimize end effects, spatial periodic 
boundary conditions are introduced. Once the functions Gfc(t) are known, the line shape 
l[E) is evaluated by means of (^. Energy will be measured in units of T whereas time 
will be expressed in units of T^^. Since we are mainly interested in the effects due to 
the presence of correlated disorder in the random system rather than in the effects of the 
different parameters on the optical absorption process, we will fix the values of V, V and 
T, focusing our attention on the defect concentration c. Thus, we have set K = 4, ^' = 10 
and T = — 1 henceafter as representative values. The width of the instrumental resolution 
was a = 0.5. The maximum integration time and the integration time step were 16 and 
8 X 10^^, respectively; larger maximum integration times or smaller time steps led to the 
same general results. The defect concentration c ranged from 0.01 up to 0.9 and for each 
lattice realization a random distribution of dimer defects was chosen. In addition, lattices 
with unpaired defects have been studied and compared with lattices containing the same 
fraction of paired defects. This enables us to separate the effects merely due to optical 
absorption in one-dimension from those which manifest the peculiarities of the correlation 
between random defects. 

In absence of defects the spectrum is a single Lorentzian line centered at E = V + 2T, 
which with our choice of parameters is E = 2. When defects are introduced in the lattice, 
a broadening of this main line is observed accompanied by a shift of its position to higher 
energies on increasing defect concentration. Keeping these general results in mind we now 
proceed to discuss the main features of the spectra obtained when defects are present in the 
system. The obtained results are shown in Figs. ^|-|^, which respectively correspond to an 
increasing defect concentration. From a close inspection of these figures several conclusions 
can be drawn. First of all, we observe that, in addition to the main absorption line located 
close to E = 2, several satellites appear in the spectra. These satellites become more intense 
as the defect concentration increases. Such lines are related to the presence of defects 
in the lattice, and their positions only depend on the particular kind of disorder (paired or 
unpaired). On increasing defect concentration, the high-energy lines increase at the expenses 
of the main line (compare Figs. [^-0). For low defect concentration (see Fig. |lD, two lines 
appear at about E = 9.0 and E = 10.2, whose positions are almost independent of c. The 
line at E = 9.0 is more intense for lattices with dimer defects whereas the line at E = 10.2 
is weaker in that case. In addition, a small shoulder appers at about E = 3.8, and it 
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becomes more intense in the case of unpaired lattices. On increasing defect concentration, 
this shoulder is clearly resolved and a new line develops at around E = 8.5 (see Fig. 
The intensity of this new line rapidly increases with c and shifts to the low-energy region 
of the spectra (see Fig. until it is placed close to E = V + 2T {E = 8 with the choosen 
parameters) . 

Firstly we discuss the shift of the main absorption line toward higher energies on increas- 
ing defect concentration. In the case of unpaired defects this shift AE depends linearly on 
defect concentration c; using a least squares fit we have found that AE ^ 0.92c in the range 
c = 0.01 up to c = 0.4. This behavior is in agreement with the average- T-matrix approxima- 
tion (ATA), which predicts that the shift increases linearly with defect concentration in the 
absence of correlations.! On the contrary, pairing of defects causes a quadratic dependence 
on the defect concentration of the form AE ^ 1.20c^. It is worth mentioning that the line 
shift AE for the paired correlated lattice is always smaller than AE for the unpaired one, for 
the same defect concentration c in the range considered. This is an important result since it 
suggests that pairing of defects reduces the effects of disorder on the Frenkel-exciton dynam- 
ics. The physical explanation of this behavior is based on the particular spatial distribution 
displayed by the defects when pairing is present. It becomes clear that the average length of 
segments which are free of defetcs is larger when correlation is present since defects always 
appear in pairs. Therefore, pairing of defects implies a smaller disruption of the lattice and 
consequently the random distribution of such defects causes less influence on the exciton 
dynamics. 

Now we focus our attention on the satellite lines present in the spectra. In order to study 
the results from a quantitative point of view, we consider the following two center problem, 
which describes the optical absorption spectrum of two isolated, coupled defect sites 

4Gi(t) = ViGi(t)+TG2(t), 
at 

t^G2it) = V2G2{t)+TG,it). 
at 

These equations of motion can be easily solved using the initial conditions 6*1(0) = 6*2(0) = 
1. Inserting the result in (0) and taking a = for the sake of simplicity, one obtains 

I{E) = I+6{E - E+) + I_6{E - E_), (4) 

where the intensity of each component is given by 

V, - ^ 



2 2 



1 + 



2T 



21 -1/2 

(5) 



V1 + V2 , ^ I fVi-V2^ 2 



and centered at 

E,^l±^±Tp^[ll^). (6) 

From the Eqs. (^ and (^) we obtain the following possible scenarios. If the two transition 
frequencies of defects are Vi=V2=V' , the intensity /_ vanishes due to the fact that only 
the dipole matrix element involving symmetric states contribute to the optical absorption 
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spectrum. In this case the spectrum exhibits a single hne centered at = V' + T = 9.0. On 
the contrary, if the transition frequencies are different (say V\ = V and V2 = V'^ the optical 
spectrum presents two components centered at Ej^ = {V + V')/2±T^1 + (V^ — V'^/AT^ = 

7 T so that E_ ~ 10.2 and E+ ~ 3.8. 

With this analytical results at hand, the quantitative interpretation of the obtained 
spectra is straightforward since each absorption line can be related to a specific kind of 
defect in the lattice. We shall consider, in first place, the case of low defect concentration 
(Fig. |l|). Besides the main line close to E = 2, characteristic satellites are clearly visible 
at about E = 10.2 (more intense in unpaired lattices) or E = 9.0 (more intense in paired 
lattices). These values agree very well will our theoretical predictions shown above: The 
line at E = 10.2 is related to centers Vi = V and V2 = V, and the number of such center 
is greater in unpaired lattices, whereas the line at E = 9.0 is related to dimer defects and 
the probability of finding such defecs in unpaired lattices is actually very small at low c. 
Analogous comments hold for the line centered at = 3.8 since it originates from the 
same defects as the line centered aX E = 10.2. On further increasing defect concentration to 
intermediate ranges (c ~ 0.5) these characteristic lines continue to dominate the high energy 
region of the absorption spectra (Fig. |^). Ultimately, in the limit c 1, the absorption 
spectrum is dominated again by a single line centered at E = V + 2T = 8, as it corresponds 
for a lattice enterely composed of defect sites. This trend is clearly observed in our spectra 
for both unpaired and paired lattices (Fig. |]). 

To seek further confirmation of the above explanations, we have investigated the intensity 
dependence on the defect concentration for the lines centered at E = 9.0 and E = 10.2. 
Results are shown in Fig. |[ The intensity of the line at E = 9.0 increases almost linearly at 
low defect concentration in the case of paired lattices, as shown in Fig. Ka). This is to be 
expected since such defects are introduced on purpose in the lattice. On the contrary, the 
intensity of this line increases quadratically with defect concentration in unpaired lattices 
since the probability of finding dimer defects in such uncorrelated lattices also increases as 
c^. Concerning the line centered at E = 10.2, its intensity increases sublinearly with defect 
concentration in the case of unpaired lattices, as shown in Fig. ^(b). This is due to the fact 
that on increasing defect concentration dimer defects appear and, consequently, the number 
of centers V\ = V and V2 = V cannot increase linearly. In the case of paired lattices 
the intensity of this line increases almost linearly since the number of centers V\ = V and 
V2 = V equals the number of dimer defects whenever the probability of pairing of such 
defects remains small at low defect concentration. Therefore, we are led to conclude that 
our interpretation of all relevant lines in the optical spectra is correct, and those spectra 
show indeed clear differences when coming from paired or unpaired disorder. 

IV. CONCLUSIONS 

In summary, we have studied the absorption spectrum corresponding to the Frenkel- 
exciton Hamiltonian for a random system in which a correlation length of the order of 
the lattice spacing is intentionally introduced, so that defects can only appear in pairs. 
By comparing the obtained spectra with that corresponding to the uncorrelated random 
system we are able to identify characteristic absorption lines which can be directly related 
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to the presence of dimer defects. Furthermore, we have obtained analytical expressions 
which explain our synthetic spectra very well, and relate microscopic system parameters, 
like transition frequencies or defect concentration, to experimental data, like position of the 
hues. Our results show that the presence of correlation effects in random systems can be 
readily determined from analysis of the absorption spectra of these samples, following the 
lines of reasoning presented in this work. 
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FIGURES 



FIG. 1. Absorption spectra for one dimensional random lattices with concentration of paired 

defects c = 0.1 (solid line) and c = 0.2 (short-dashed line). Results arc compared to lattices with 
the same concentration of unpaired defects c = 0.1 (long-dashed line) and c = 0.2 (dot-dashed 
line). All spectra have the same area. Note the magnification factor of the high-energy part of the 
spectra. 

FIG. 2. Absorption spectra for one dimensional random lattices with concentration of paired 
defects c = 0.3 (solid line) and c = 0.4 (short-dashed line). Results are compared to lattices with 
the same concentration of unpaired defects c = 0.3 (long-dashed line) and c = 0.4 (dot-dashed 
line) . 

FIG. 3. Absorption spectra for one dimensional random lattices with concentration of paired 
defects c = 0.7 (solid line) and c = 0.9 (short-dashed line). Results are compared to lattices with 
the same concentration of unpaired defects c = 0.7 (long-dashed line) and c = 0.9 (dot-dashed 
line) . 

FIG. 4. Intensity of the absorption lines centered at a) E = 9.0 and h) E = 10.2 as a function 
of dimer defect concentration in random one-dimensional lattices with paired defects (full circles) 
and unpaired defects (full triangles). Solid lines indicate the best least squares fit. 
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